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using partition of qudit sytem into a set of qubit system 
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Gerjuoy [Phys. Rev. A 67, 052308 (2003)] has derived a closed- form lower bound for the entanglement 
of foiTnation of a mixed qubit-qudit system (qudit system has d levels with d > 3). In this paper, inspired by 
Gerjuoy 's method, we propose a scheme that partitions a qubit-qudit system into d(d— 1)/2 qubit-qubit systems, 
which can be treated by all known methods pertinent to qubit-qubit system. The method is demonstrated by a 
qubit-qudit system (The levels of qudit are d = 3 and d = 5, respectively). 
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I. INTRODUCTION 



Entanglement of quantum systems is an important physi- 
cal resource to realize quantum information tasks and quan- 
tum computation The quantitative measure of entan- 
glement is one of the main research areas in quantum infor- 
mation theory and quantum computation [2], which have at- 
tracted much attention of many researchers. In this sense, 
many useful measures were developed, such as: concurrence 
ol- Uin . entanglement of formation (EOF) |12Ul7ll . geomet- 
ric measure [ 18-2Tl, entanglement witness f22,'23l, quantum 
discord [24, 251 , three-tangle [26J, etc. These entanglement 
measurements are usually defined first for pure states and then 
are extended to mixed states via the convex roof construction. 
Because it requires complicated optimization procedure, gen- 
erally speaking, computing an entanglement measurement for 
a given mixed quantum state is difficult. To the best of our 
knowledge, only a few analytic expressions of some entan- 
glement measurements for some special quantum systems, for 
example, the EOF and geometric measure of qubit-qubit states 
and isotropic states, are 

obtained llllll mi H- 

some numer- 
ical algorithms for computing some entanglement measures 
were also developed [201 . 

The concuiTence and the EOF among all entanglement 
measurements play an essential role due to some of other en- 
tanglement measurements can be expressed by concuiTence 
and the method employed to derive analytic expressions of 
the EOF can be used to derive analytic expressions of other 
entanglement measurements. Wootters have obtained an el- 
egant formula for qubit-qubit system Wei and Goldbart 
have also derived an analytic expression of the geometric mea- 
sure for two-qubit mixed states [18]. Due to the fact of the 
concurrence, as defined originally, is only suitable for qubit- 
qubit systems and the optimization process to get analytic ex- 
pressions for the EOF of a general entangled state in higher- 
dimensional space to be complicated, several schemes were 
proposed to find the lower bounds of the concurrence and the 



EOF of general entanglement mixed states Using 
the Schmidt decomposition theorem [H^, Geijuoy derived the 
lower bounds of the concuiTence and the EOF of any qubit- 
qudit system, and easily obtained Wootters' formula. The 
method employs a set of (2 x d) x (2 x d) matrixes S*-*, that 



(2) LL ^ LL 

are constructed via ay €5 <Ty"^ where ay'^' (cry"'') is the usual 
ay Pauli matrixes for the qubit (qudit) case. In this sense, in 
the present paper we propose a scheme that partitions a qubit- 
qudit system into d{d — l)/2 qubit-qubit systems, which can 
be treated by all known methods pertinent to qubit-qubit sys- 
tem. 

This paper is organized as follows. In the next section, 
we analyze Gerjuoy's method and propose our scheme, which 
simplify Gerjuoy's procedure. In the section|II] two examples, 
qubit-qutrit (d = 3) and qubit-qudit ( with d = 5) systems, are 
given to illustrate our scheme. The summary and discussion 
are given in section HV] 



(2) rJd). 



II. ANALYSES OF GEFIUOY'S SCHEME AND 
PARTITION OF QUBIT-QUDIT SYSTEM INTO 
QUBIT-QUBIT SYSTEM 

To find the lower bound of entanglement of formation 
(EOF) of the qubit-qudit system, E. Gerjuoy |6] first defined 
d{d — l)/2 symmetric square matrices S'-' , < z < d — 2 and 
j > i, whose elements all are zero, except for 



^i+d,j 



1. 



Second, he defined 



Cij{p) = max{0, A*/ - X'^ - A^f - X''^), 



(la) 
(lb) 

(2) 



where the A'^ , ordered decreasingly, are the square roots of 
the four largest eigenvalues of the matrix pS*-'p*S'-', p is the 
density matrix of the qubit-qudit system and p* is its conju- 
gate. Thirdly, he denoted the lower bound of the concurrence 
C(p) of the qubit-qudit system by Cdb{p), 
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j>i i=0 



1/2 



< C{p). 



(3) 



2 



The desired lower bound on the qubit-qudit EOF is e [Cdb {p)]- 
For a qubit-qutrit mixed state system, there are only three 
S'-' expressed as Sx,Sy and Sz, 
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(4) 



(5) 



(6) 



respectively. For the qubit-qubit mixed state, there is only one 
5*-' denoted as 



S = 



-1 
10 
10 
-10 



(7) 



Note that the Eq. (|7]l is a matrix constructed via cr^ (g) cr^, with 
A and B denoting the indexes of qubit A and qubit B, respec- 
tively. Considering a mixed qubit-qubit quantum system with 
levels |0) and |1) one can construct the PauU operator ay for 
the subsystems A and B given by 



of cr^oi ® '^yij °'^ly ^ minus. It doesn't affect the 

eigenvalues of pS'-'p*S*-' if we replace S'-' by cr^oi ® '^y^ij- 

It is now clear that Gerjuoy's approach, in fact, is to treat 
the qubit-qudit mixed system as a set of (i!/(2!(d — 2)!) qubit- 
qubit system. Therefore, for such a system, we can take any 
two levels of qudit subsystem to combine with the qubit sub- 
system as a qubit-qubit. Then, all methods to solving the 
qubit-qubit problem can be used. In our present case, we only 
need to use matrix S instead of S*-* (or S^; , Sy and Sz) to calcu- 
late the concurrence C'-' (C^; , Cy and Cz) for those qubit-qubit 
subsystem. This will greatly simplify the calculation. In the 
next section, we shall demonstrate our proposal by a concrete 
example. 



III. ILLUSTRATIVE EXAMPLES 

To illustrate our method let us consider two atoms (A and 
B), each of them interacting resonantly with a single quantized 
mode of a cavity field (system C) in a Fock state. This phys- 
ical situation is described by the two-atom Tavis-Cummings 
(TC) Hamiltonian: H — hg[{aA + crB)ap + (cr^ + (T^g)ac], 
where CTj and crj are the Pauli ladder operators for the jth 
atom, a(at) is the annihilation (creation) operator for pho- 
tons in cavity C, and g is the coupling constant. We as- 
sume the system is initially in the state |?/'(0)) — (a|0^0^) + 
j3\lj^lg))\n^). Since the TC Hamiltonian preserves the total 
number of excitations, the cavity mode will evolve within a 
five-dimensional Hilbert space spanned by { |(n — 2)^ ) , | (n — 



l)c), 



|(n + l)c), |(n + 2)c)}forn > 2. When n = 0, 1 



j(|lx)(Ox|-|0,)(l,|), {x^A,B) 



(8) 



the dimension will be 3 and 4, respectively. On the other 
hand, the atomic system will evolve within the subspace 
{|0,0J, |+), |1,1 J} with 1+) = (|1_,0J + |0,1 J)/^/2 
independently of n. By solving the Schrodinger equation, the 
system at time t is described by the state 

\m) = Ci(i)|0,0j|(n + 2)^>+C2(i)|+)|(n + l)^) 
+C3W|l^l,|nJ+C4(t)|0^0j|n^) 
+C5W|+)|(n-l)J+C6(t)|l,l,|(n-2)^), 

(10) 



Then, if we consider a mixed qubit-qutrit system composed of 
subsystem A with two levels |Oyi) and |lyi) and a subsystem 
B with three levels |0b), |ls), |2b), the Pauli operator ay can 
be now expressed by Dirac notation as: 



i{\ix){iA-%){3x\). {x = A,B), 



(9) 



where = {0, 1} for a; = A and = {0,1}, {0,2} 

and {1, 2}, respectively, for x = i?. It is easy to test the matrix 
forms of CT^oi^cr^oi' '^ym®'^y02 ™d crym^dy 12 in basics 
|0a,0b), |o'^, 1b),'|0^, 2b), |U',Os>, 1b), li^, 2b) are 
Sx, — S„ and — Sz 



where the probability amplitudes are 



ciW - o„ , o ^[1 - cos(V2(2n + 3)5t)], 



2n + 3 



C2{t) = 



i|3^/n + 1 
V2n + 3 



qubit-qudit mixed state system the matrix forms of cr^oi 



{i < j) in basics |0a,0b), |0a,1b), -, |0a, - 
|1^,1b), |lA,(d-l)s) equals -S'^ 



1)b), |1a,0b;, |iA,-LB; 
It is notable that the difference between the matrix expression 



C4(i) = a < 1 — 



2n - 1 



(11) 



sin(V2(2n + 3)gt), (12) 



A,^B/, |-^A,UB/, iB/, |J^A,^B/ lie ^ n+l 1 

respectively. Generally speaking, for a ^{t) = 13 ll - -[1 - cos(V2(2n + 'i)gt)] > , (13) 

le svstem the matrix forms of it^„, (X) 2n + 3 J 



[l-cos(V2(2n-l)5t)]L (14) 



3 



C6(t) 



C5 (0 = sin(V2(2n-l)5i), (15) 

V2n — 1 



y\Jn{n- 1) 
2n- 1 ' 



3(V2(2n-l).gi)]. (16) 



Now, we take trace of density operator p = \ip{t)){il]{t)\ 
over atom B resulting in the reduced density operator of the 
qubit-qudit system pAC- 




0.0 0.2 0.4 0.6 



A. Qubit-Qutrit case 

When 71 = 0, atom A and cavity C compose a qubit- 
qutrit system. As described in the above section, we delete 
terms not containing ji^, Oc)Oa, Oc|,|«A, Oc)Oa, lc|, 
|m,1c>0a,0c| and \iA,lc){jAAc\ {hi = 0,1) in pAc to 
form 

Pac = c3c;|1a,0c)(1a,0c|+c4c:|0a,0c>(0a,0c| 
+ ^C2C^(|0a, 1c)(0a, Icl + lU, 1c)(1a, Icl) 

+ -^(c2c:|1a,1c)(0a,0c|+c^c4|0a,0c)(U,1c|) 

+ -^(c3c;|U, Oc)(Oa, Icl + c^calO^, 1c)(1a, Oc|), 

(17) 



FIG. 1: (Color online) The evolution of entanglement of foima- 
tion between the atom A and the cavity mode C for the initial state 
|V(0)) = (a|0^0^) +/3|l^l^))|0c) withQ = p = 1/V2. The 
red solid line corresponds to the present expression of Cac and the 
blue dashed line to the lower bound of Cac given by Ref. The 
dimensionless time r — \/6gt/{2iT). 



{0,0,|clc3|,|clc3|}, therefore, Cy 



0. Conse- 



quently, the lower bound of concurrence C{p) of this qubit- 

qutrit system is Cac 
The lower 

EfiCAc) 



'2[|ciC2r + (|c2C4|-|C2C3|)2]. 

The lower bound of EOF for the qubit-qutrit system is 
h{{l + yr~~C^)/2), where h{x) = 



-a;log2 



(1 — x) log(l — x) is double entropy func- 



The matrix form of which is 

/ C4C4 



P°a'c 



C4cl/V2\ 

C2C^/2 C2cyV2 

C3C*JV2 C34 

\c2cl/V2 C2C^/2 / 



Similarly we obtain matrices p^^ and p]^(j respectively 



(18) 



p'a'c 



I C4CI C4CI 

cic| cic| 

C34 

V 



(19) 



12 
Pac 



( C2C^/2 




V 







\ 



CiCi 

C2C^/^/2 c2C^/2 
/ 



(20) 



The matrices p^^^ and are X form [28]. The correspond- 
ing concurrences can be read out 



Ca; — Coi — V2 I 



C2C4I 



IC2C3I 



(21) 



tion. The evolution of Eac with the dimensionless time 
T = \/&gt/{2T:) are plotted in Fig.l. In order to make a 
comparison, Eac also was computed according to the lower 
bound of Cac given by Ref. Ist]. Though two lines in Fig.l 
do not coincide, but their behaviors of evolution with time are 
the same. 



B. Qubit-Qudit case 

Now, considering n = 2, the atom A and cavity C compose 
a qubit-qudit {d — 5) system. According to our scheme, we 
can take any two levels of cavity mode C to form ten qubit- 
qubit system with two levels of atom A and find: 



p"2c 
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C5C5 
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(23) 
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C3CI 


C3C3 



(24) 



= C12 = V2 



IC1C2I 



(22) 



Unfortunately, The matrix p^^ is not the X form. 
The square roots of four eigenvalues of p^^jSp^^S are 
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(31) 
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(27) 
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(32) 
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— 6— 



C4e, \ 






(29) 



(30) 



The eigenvalues of p^^Sp^^S are all equal 

to that of pAcS'^ p*Ac^'^ (« = 0,...3, j > 
i), The non-zero square roots of them are 

{V2\c5Ce\}, {\c4Ce\, \c4Ce\}, {V^\c2Ce\},{\ciC6\,ciCe\}, 

{y2|c3C5|, y2|c4C5|}, {V2|ciC5|}, {V2|C2C3, V2\c2Ci\}, 

{IC1C3I, IC1C3I}, {-\/2|ciC2|}, respectively. Correspondingly, 

Coi = V2\c5Ce\,Co3 = V2\c2Ce\,Ci2 = V2||c3C5| - 

|c4C5||,Ci4 = \/2|ciC5|, C23 = y2||c2C3| - |c2C4||,C34 = 

■\/2|ciC2|, C02 = Co4 — Ci3 — C24 = 0. The lower bound of 
concurrence C(p) of the qubit-qudit is 



Cb5 = V2J (IC2C3I - |C2C4|)2 + (IC3C5I - |C4C5|)2 + cfc| + cfc| + c|c§ + ci 



,2„2 
5^6- 



(33) 



r 



We plot the evolution of Eac with the dimensionless time 
T = ^/Ugt/ (Gtt) in Fig. 2. i^yic is also computed according 
to the lower bound of Cac given by Ref. |9] and plotted in the 
same figure for comparison. Though two lines in Fig. 2 have 
different trends in some interval of r, their global behaviors 
of evolution with time are basically the same. 



IV. SUMMARY 

We have analyzed Gerjuoy's approach on calculating the 
lower bound on entanglement of formation for qubit-qudit 
system and we found that his method, in fact, is to treat qubit- 
qudit system as a set of qubit-qubit system. Therefore, we 
proposed a simple scheme to solve qubit-qudit problem. The 
scheme consists of three steps: (1) partition the qudit system 
into a set of qubit system; (2) compose the original qubit and 
partitioned qubit into a set of qubit-qubit systems and treat 
them by all methods suitable to qubit-qubit system. Find the 



measurements you want for every qubit-qubit system; (3) ob- 
tain the measurement of whole qubit-qudit system. For the 
case discussed in the present paper, we calculated the concur- 
rences for every qubit-qubit system and the lower bound of 
the concurrence of the qubit-qutrit or qubit-qudit system. Our 
method has the advantage of avoiding finding many matrices 
S*-' and only using one matrix S — cjy ® <jy. This method 
greatly simplified the calculation about the measurement of 
qubit-qudit system. We hope this method can be extended to 
treat other problems of the qubit-qudit system. 
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FIG. 2: (Color online) The evolution of entanglement of forma- 
tion between the atom A and the cavity mode C for the initial state 
|V(0)) = (q10^0^)+/3|1^1^))|2^) with a = /3 = 1/V2. The 
red solid line corresponds to the present expression of Cac and the 
blue dashed line to the lower bound of Cac given by Ref. Q]. The 
dimensionless time r — \/lAgt/{6Tv). 
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